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Abstract 



A microscopic calculation of near-threshold //-meson production in the reac- 
tion NN — ► NNn is presented. It is assumed that the r\ meson is produced 
via direct emission and via elementary rescattering processes MN — ► r/N of 
various mesons M = ir, p, etc. As a novel feature the amplitudes for the 
latter production mechanism are taken from a multi-channel meson-exchange 
model of the ttN system developed by the Julich group which contains ex- 
plicitly the channels ttN, pN, rjN, aN, and 7rA. Furthermore, effects of the 
NN interaction in the final as well as in the initial state are taken into ac- 
count microscopically. Our results are compared with recent data from the 
COSY and CELSIUS accelerator facilities. Reasonable agreement with avail- 
able near-threshold cross section data for the reactions pp — > ppw, pn — > pun, 
and pn — ► dn is found. 
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I. INTRODUCTION 



Meson production in nucleon-nucleon (NN) and nucleon-nucleus collisions is a rather attrac- 
tive research field in hadron physics at intermediate energies. The revived interest in such 
reactions is closely connected to recent successful developments in the accelerator technology 
which enabled experimentalists to perform rather accurate measurements on the production 
of 7r, n, n' and other mesons in nucleon-nucleon collisions in the near-threshold region [1-3]. 
These data in combination with the theoretical analysis open the unique possibility to study 
the mechanism of the production of different mesons as well as to produce information on 
resonance properties, baryon-baryon and meson-baryon interactions, short range properties 
of the NN interaction etc.. 

Among those reactions the production of the n meson is certainly rather interesting. E.g., 
the rj is the only other non-strange pseudo Goldstone boson existing and thus a close though 
heavier relative of the pion. Also it is expected that n production near threshold is closely 
linked with the properties of the Su A r *(1535) resonance. Finally, and may be most im- 
portantly there is already a wealth of rather accurate data on n meson production near 
threshold providing a sensible testing ground for model calculations. Apart from total cross 
sections for the reaction pp — > ppn [4-9], there are also some data on the total cross sections 
of the reactions pn — > pni] [10] and pn — > drj [11,12]. Furthermore, there are data on differ- 
ential cross sections and invariant-mass spectra [9,13-15] and even anlyzing powers for the 
reaction pp — > ppn [16]. 

The present work is devoted to the study of the production of rj mesons in NN collisions. 
The goal is a combined theoretical analysis of all measured channels of the reaction NN — > 
NNrj, namely pp — > ppi], pn — > pni] and pn — > dn, by taking into account consistently the 
interaction between the nucleons in the final as well as in the initial state and by utilizing 
a microscopic model of meson-nucleon (MN) scattering [17,18] for the description of the 
77-meson production process. We focus on the description of 77-meson production in the 
near-threshold region, i.e. for excess energies Q up to around 50 MeV. Therefore, we will 
restrict our investigation to S-wave contributions. Recent data on differential cross sections 
of the reaction pp — > ppn [14] suggest that in this energy region the production of n mesons 
still takes place mainly in S-waves. 

In spite of the abundance of theoretical investigations on ?i-meson production [19-32] the 
leading production mechanism is still not identified. There is general consent in the literature 
that 77- meson production is dominated by re-scattering processes as depicted in Fig. la, 
with a variety of intermediate meson exchanges. However, it is not clear which one of the 
possible meson exchanges plays the most important role. For instance, in Refs. [21,22,25,26] 
it is suggested that the dominant contributions should come from the p meson, whereas in 
[27,30] it is found to be of minor importance compared to other contributions and in [24] that 
particular contribution is not even considered. Recently Nakayama et al. [32] compared the 
two scenarios of a n- and a p-exchange dominance and concluded that both these scenarios 
can describe the cross section data for the various NN — > NNrj channels equally well. 

Apart from the choice of the exchanged mesons there are also differences in the construction 
of the MN — > nN transition amplitude. Many of the models [19-22,25,27,32] are based 
on the calculation of tree level re-scattering diagrams where the MN — > nN transition am- 
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plitude is simply parametrized by the Sn iV*(1535) resonance. Background contributions, 
e.g., through t-channel meson exchanges are either omitted or simply added on top of the 
resonance contribution. In such approaches the relative phases between different meson 
exchanges have to be fixed by hand, which naturally introduces large uncertainties. As a 
consequence quite contradictory prescriptions have been adopted in the literature. E.g., in 
Ref. [21] it was tried to minimize interferences between the various production mechanisms. 
Assuming that the 77's are produced only via N* excitation led to 7r and p exchange contri- 
butions that are orthogonal to each other and, therefore, don't interfere at all. On the other 
hand, in a recent investigation by Faldt and Wilkin [26] interferences are sort of maximized. 
Those authors argue that only a strong interference of the p exchange amplitude with the 
it exchange contribution allows to achieve a consistent description of the experimentally 
observed cross sections for the reactions pp — > pprj and pn — > pni]. 

In the present investigation we employ MN — > r]N amplitudes generated from a microscopic 
(meson-exchange) model of the ttN system. It is a coupled-channels model that contains 
5 channels, namely irN, rjN, pN, aN, and 71 A [18]. The interactions in and between the 
various channels are derived in the meson-exchange picture starting out from effective chiral 
Lagrangians. The model includes the iV*(1535) resonance as essential contribution but in 
addition also various (t-channel) meson and (-u-channel) baryon exchange diagrams. The 
parameters of the model are fixed by requiring a simultaneous description of the ttN phase 
shifts and inelasticities as well as of the transition cross sections for the reactions ttN — > rjN 
and irN — > pN over an energy range that extends well beyond the r/N threshold. The reac- 
tion amplitudes are obtained from solving a coupled-channels relativistic scattering equation 
of Lippmann-Schwinger type. Clearly, in such a model not only the nN — > rjN amplitude is 
determined by empirical data, but also the transition amplitudes involving heavier mesons 
are to a large extend constrained by the phase shifts and inelasticity parameters of ixN 
scattering. This is also true for the relative phases between the various amplitudes. 

Another uncertainty in investigations of the reaction NN — > NNrj is caused by the treatment 
of the NN FSI and ISI effects. In many papers concerned with meson production FSI effects 
are taken into account approximately by applying methods which are rather questionable if 
one wants to obtain absolute predictions for the production cross sections (see, e.g., Refs. 
[19-21,23,27] for the case of 77-production), as was pointed out recently in Refs. [33,34]. The 
situation with the ISI is even less satisfying. Indeed, its effect is simply omitted in most 
studies. 

In our investigation we take the FSI and ISI in the NN system fully into account. In 
particular, we employ an NN model that reproduces the relevant NN phase shifts reasonably 
well up to energies around the r\ production threshold and, therefore, allows a consistent 
description of the FSI and ISI effects. 

Let us mention here that, at the present stage, the Jiilich MN model does not include 
the ujN channel. Thus, we cannot take into account the cu-rescattering contribution to the 
reaction NN — > NNrj. Anyhow, we want to emphasize that most model studies indicate 
that its contribution should be small [20,21,26,32]. 

The paper is structured in the following way: In Sect. II we describe the ingredients that en- 
ter into our model calculation of 77 production in nucleon-nucleon collisions. In particular, we 
present some details about the meson-nucleon amplitudes that are utilized for constructing 
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the elementary n meson production operator and we specify the relevant vertex parameters. 
We also provide a brief description of the NN model that is used for the treatment of the 
ISI and FSI, and we supply its prediction for those energies and partial waves in the initial 
NN state, relevant for n production near threshold. Our results are presented and discussed 
in Sect. III. We show our predictions for various charge channels, i.e. pp — » ppi], pn pur], 
and pn — > dn, and compare them with available data on total cross sections. We also study 
the sensitivity of our results to variations in the elementary (pN — > nN) production ampli- 
tude and to differences in the NN interaction in the final state. Finally, we analyse the role 
played by the various production mechanisms included in our model. A summary and some 
concluding remarks are given in Sect. IV. 

II. MODEL CALCULATION OF THE REACTION NN -» NNrj 

We study r^-meson production in distorted wave Born approximation (DWBA). This means 
that the transition amplitude from the initial to the final state, M 2 ^ 3 , is given by the 
expression 

M^ 3 = (Xf I (1 + TP N I Go)A 2 ^ 3 (l + G TlM) \ Xi ). (1) 

where is the (elementary) production operator and T^j and T^f/ are the NN reaction 
amplitudes in the initial and final states. A graphical representation of Eq. (1) is given in 
Fig. 2. Note that here we have tacitly assumed that the final-state effects are dominated 
by the interaction in the NN system and that the interaction in the nN system can be 
neglected. There is, however, experimental evidence [6,11,35] that the nN interaction has a 
significant impact on the energy dependence of the cross section. Under such circumstances 
it would be, in principle, desirable to include the nN FSI consistently and to solve 3-body 
integral equation of Faddeev type for the most difficult case when all particles are in the 
continuum (see, e.g. Ref. [36] and references therein). A few efforts along this line can be 
already found in the literature [37,38], though only for the quasi two-body reaction np — > dn. 
On the other hand, the experiments also indicate that the nN interaction might have an 
influence only for very small excess energies, i.e. up to around 10-15 MeV above the n- 
threshold. Thus, a sensible investigation of the 77-meson production mechanisms should be 
still feasible without taking into account the nN interaction, if one focusses on the energy 
range 10 < Q < 50 MeV, say. In any case it is quite reasonable to consider the much more 
complicated 3-body effects only after a detailed quantitative understanding of NN ISI and 
FSI effects as well as of the production process has been achieved. Therefore, in this work 
we disregard possible effects from the nN FSI. 

The production amplitude A2-»3 in our model consists of the re-scattering diagrams with 
7r, p, n, a meson exchanges and the direct n production. Corresponding diagrams included 
in our work are shown in Fig. 1. 

One of the principal novelties of our investigation is the utilization of a realistic microscopic 
model for the elementary reaction amplitudes MN — > nN, namely a coupled-channels model 
for ttN scattering that has been developed recently by the Jiilich group [18,39]. The im- 
portant aspect connected with the use of this model is that the off-shell properties of those 
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amplitudes and also the relative phases between different meson contributions are now fixed 
and (to a certain degree) constrained by the data on ttN scattering. Another merit of our 
model study of r\ production in NN collisions consists in the full and consistent treatment 
of effects from the NN FSI as well as ISI. 

As mentioned already we restrict our investigation to the S-wave contributions. Recently 
measured angular distributions in the reaction pp — > ppr\ exhibit an isotropic structure at 
the excess energies Q — 15 MeV as well as at Q = 41 MeV [14], thus suggesting that even 
up to Q ~ 50 MeV the production of 77 mesons might take place predominantly in S-waves. 
We want to point out, however, that the invariant-mass distributions reported in the same 
paper cannot be described by S-waves plus pp FSI alone [14]. This could be an evidence 
that at least at the higher energy P- waves already play a role [40] . 

Restricting ourselves to S-waves we have only two amplitudes corresponding to the possible 
final states ^^SqS (for isospin 1=1) and 3 SiS (for isospin 1=0). The corresponding initial NN 
states are 3 Po an d 1 Pi, respectively. (We use the standard nomenclature where capital letters 
denote the NN partial waves and the small letter labels the orbital angular momentum of 
the 77 with respect to the NN system.) 

A. Model of the meson-nucleon (MN) interaction 

In this subsection we describe the structure and results of the MN model - which is one of 
the main ingredients of our calculation. This model was developed in Refs. [17,18,39]. The 
ambitious aim of the model is a simultaneous description of the nN phase shifts and inelas- 
ticities as well as of the transition cross sections for the reactions ttN — > i]N and itN — > pN . 
The reaction amplitudes are obtained from solving a relativistic coupled-channels scatter- 
ing equation of Lippmann-Schwinger type. The model includes 5 channels, namely ttN, 
r/N, pN, aN, and irA. All the MN interactions and transition potentials are constructed 
in time-ordered perturbation theory (TOPT) on the basis of effective chiral Lagrangians. 
Most of the coupling constants at the meson-baryon-baryon and meson-meson-meson ver- 
tices included in the model are taken from other sources [18]. The cutoff masses in the 
vertex form factors are the only free parameters which were determined by a fit to the ixN 
partial-wave amplitudes. The potentials consist of t-channel meson exchanges, s-channel 
resonance graphs and -u-channel baryon exchanges (cf. Figs. 2-5 in Ref. [18]). 

The model provides a satisfactory overall description of all ixN partial amplitudes with 
total angular momentum J < 3/2 [18]. Specifically, in the Su partial wave, which is the 
relevant one for the reaction NN — > NNrj close to threshold, the agreement with the 
experimental information in a region of about 60 MeV above the 77 threshold (the r]N 
threshold is at E cm 1486 MeV) is quite good, as demonstrated in Fig. 3 by the dashed 
lines. Therefore this model is a reasonable starting point for our calculation. The off- 
shell transition amplitudes of this model are utilized for the meson rescattering amplitudes 
MN — > r/N that enter into the evaluation of the production operator A 2 ^3 (cf. the filled 
circle in Fig. la). 

It should be said, however, that the data on ttN scattering (plus the irN — > r]N and 
iiN — > pN transition cross sections) do not fully determine all the transition potentials 
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of the coupled-channels model. E.g., in the original model [18] no direct pN — > i]N tran- 
sition potential at all was included, because the 7rN results are not too sensitive to such 
contributions and, accordingly, corresponding parameters could not be fixed unambigously. 
Possible contributions to pN — > r/N can arise from t-channel meson-exchange (e.g., p, b 1: 
etc.), from N exchange and, in particular, from the s-channel pole diagram involving the 
Su A r *(1535) resonance, as shown in Fig. 4. 

Note that the original 7rN model [18] still generates a pN — > r]N transition amplitude - even 
in the absense of a direct pN — > rjN transition potential. Such an amplitude arises naturally 
in a coupled-channels model, e.g. via higher order transitions of the kind pN — > ttN — > r)N. 

In view of the importance of the pN — > 7]N transition amplitude in many of the earlier 
investigations of r\ production, cf. the Introduction, and specifically in view of the crucial 
role it plays in the analysis of Faldt and Wilkin [26] we want to explore how well this 
amplitude is constrained within our coupled-channels model and to which extend it can be 
varied. For this purpose we have created a variant of the MN model where we now explicitly 
include a coupling of the pN channel to the Su iV*(1535) resonance, cf. the diagrams in 
Fig. 4. Evidently, we need to make sure that the description of the nN data in the region 
not far from the r\ threshold remains basically unchanged when those diagrams are added 
to the model. The bare coupling constants g w NN*, Q^nn*, and g P NN* have been varied to 
explore the variation in the amplitude for pN — > rjN while staying as close as possible to 
the experimental phase shift and inelasticity of the Su ttN partial wave (cf. Fig. 3) and the 
nN — > r]N transition cross section produced by the original model (cf. Fig. 5) in a region 
of about 60 MeV around the r/N threshold. 

All other nN partial waves remain unchanged since the parameters of the t- and u-channel 
exchange potentials were not varied. 

Eventually we settled on a small and negative value for the bare coupling constant g P NN* 
since it turned out that only choosing it to be negative allows to achieve a significant influence 
of those additional diagrams on the pN — > rjN amplitude and, in turn, on the predictions 
for i] production in NN — > NNrj reactions. In Sec. Ill we present results for this extended 
model together with those based on the original MN model. 

As is evident from Fig. 5 the present MN model of the Jiilich group overestimates the 
n~p — > r/n cross section by roughly 15%. In order to estimate the impact of this shortcoming 
on the results for NN — > NNt] we performed some exploratory calculations based on a 
reduced ttN — > rjN amplitude. Specifically, we weakened the ttN — > r/N transition potential 
of the MN model [18] phenomenologically by decreasing the coupling of the a meson at 
the a TT7] vertex [42] so that the transition cross section is reproduced. (Note that the irN 
phase shifts are no longer described in this case!) It turned out that a decrease in u^-p^n 
by 15% implies a reduction in the predictions for u nn—>nnt] by roughly the same amount. 

B. Vertex parameters 

Besides the MN — > r]N amplitude the production amplitude ^2^3 also contains the meson- 
nucleon-nucleon (MNN) vertices from where the rescattered mesons are emitted, cf. Fig. 
la. In this subsection we want to provide the parameters involved in those vertices, i.e. 
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the coupling constants and the cutoff masses in the form factors. As far as the coupling 
constants are concerned all of them are taken over from the full Bonn NN model [43] . The 
only exception is the one of the 77 meson - because this meson is not included in the full Bonn 
model. Here we take the value which follows from SU(6) symmetry, i.e. gi NN /{^) = 1.8, 
which is also used in the MN model [18] for the r/NN vertex in the nucleon-exchange 
diagrams. 

The vertex form factors are likewise taken over from the full Bonn NN model or (in case of 
the f) meson) from the 71N model [18]. In those models the form factors are assumed to be 
of conventional monopole form, i.e. F(q M ) = (A 2 MNN - m 2 M )/(A 2 MNN + q^), where A M nn 
is the cutoff parameter and qu and m M are the 3-momentum and mass of the exchanged 
meson, respectively. The values of the employed vertex parameters are summarized in Table 
I. 

Since the elementary MN — > rjN amplitude and also the vertex parameters are fixed from 
earlier investigations we do not have any free and/or adjustable parameters in the produc- 
tion amplitude A 2 ^ 3 . Thus, our results for the reaction NN — > NNrj are genuine model 
predictions. 

C. NN interaction in the initial state 

The laboratory energy corresponding to the ^-production threshold is = 1250 MeV. 
Thus, the NN interaction in the initial state takes place at rather high energies. At such 
energies the effects of the ISI are characterized by the following features: 

• The ISI has practically no influence on the energy dependence of the 77 production 
cross section because the variation of the NN interaction in the initial state within 
the energy interval determining the threshold region is negligible [44]. 

• The NN scattering in the energy region in question is already strongly inelastic because 
multiple-pion production channels are open. In Ref. [33] it is shown that in such a 
case the ISI leads to a significant reduction of the cross section. 

Indeed, the energy relevant for the ISI is so high that the well known realistic NN models 
such as the Bonn and the Paris potentials [43,45] can not be applied anymore. In case 
of the Paris or the Bonn one-boson-exchange (OBE) potentials no inelastic channels were 
included and therefore they are not valid for energies above n production threshold. In the 
full Bonn model such inelastic channels are, in principle, build in via couplings to the NA 
and AA channels. However, since meson retardation effects are retained as well three body 
singularities occur for energies above n production threshold and hence the calculation is 
technically much more difficult [46]. Therefore, in the following we use an alternative NN 
model, CCF [47], which contains the same dynamics as the full Bonn model, specifically 
the coupling to iVA and AA channels, but where meson retardation effects are removed 
by the so-called folded diagrams technique. In Figs. 6 and 7 we present the results of the 
CCF model [47] for the NN phase shifts and inelasticities in the 3 Po and x Pi partial waves. 
Since we restrict ourselves to final states composed purely out of s-waves, these are the 
relevant partial waves for the NN interaction in the inital state. The agreement with the 
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experimental data is quite good, specifically in view of the fact that the model was only fitted 
to the NN phase shifts below the pion production threshold, i.e. for Ti a b < 300 MeV. Note, 
that in this model the inelasticities are introduced only through the coupling to channels 
involving the A isobar. Therefore, in the x Pi channel, where only the simultaneous exitation 
of two A isobars is allowed, the inelasticity becomes noticable only after 600 MeV kinetic 
energy, i.e. somewhat later than indicated by the empirical data. However, in the energy 
range relevant for r) production the description is reasonable. Finally, let us mention that 
the width of the A isobar is taken into account phenomenologically by using a complex A 
mass in the propagator. Specifically, we employ a parametrization of the width given by 
Kloet and Tjon in Ref. [49] which is energy- as well as momentum-dependent. 

The evaluation of the ISI effects amounts to performing a loop integration, cf. Eq. (1) and 
Fig. 2, which gives rise to contributions from the unitarity cut and from the principal value 
(PV) integral. Since in the initial state the colliding nucleons have rather large relative 
momenta it was argued in Ref. [33] that the contribution of the PV integral should be 
suppressed as compared to the piece resulting from the on-shell (unitarity) cut. Under this 
assumption the ISI effect on the cross section reduces to a simple multiplicative factor |A| 2 
which can be expressed in terms of on-shell NN information only, namely phase shifts (5) 
and inelasticities (rj): 

\\\ 2 = |l -i^P p T NN \ 2 = hi + r] cos(25) + ir] sin(25)\ 2 . (2) 
2 4 

2 r> e 2iS — 1 

Here, the T-matrix is defined as T = and p (E„) is the cms momentum 

nE p 2ip 

(energy) of the nucleons. 

In our actual calculations it turns out that there is a strong cancellation between the Born 
term and the piece coming from the unitarity cut, i.e. between the terms corresponding to 
the "1" and to ll r]cos{28f in Eq. (2). Then the ISI factor A in Eq. (2) is rather small. But 
in such a case, the contribution from the PV integral, which is certainly small and which 
was neglected in Eq. (2), becomes important as well and has an influence on the ISI effects. 
We will come back to this point in the next section. 

Due to 3-body singularities appearing in the loop with the ISI the evaluation of the PV 
integral is very involved [50]. In the present work we avoid the technically rather tedious ex- 
plicit treatment of this 3-body singularity. Rather we follow the approach which was already 
used in the work of Batinic et al. [24] and suppress the 3-body singularity by putting the 
nucleon energies appearing in the meson propagator in the intermediate state on the energy 
shell. We want to emphasize that this approximation still provides the correct value for the 
contribution from the two body NN unitarity cut in the ISI loop. Only the contribution 
from the PV integral is influenced by the approximation discussed above. 

The effects of the ISI on the results for rj production in NN collisions will be discussed in 
detail in section III. Here we want to emphasize only that the PV integral in the ISI loop 
introduces an interesting feature. It breaks the universality of the ISI as it is suggested by 
the prescription given in Eq. (2), i.e. it implies that the reduction caused by the ISI will 
differ for the various contributions to 77-meson production and, thus, depends also on the 
production mechanism. The approximative treatment of ISI according to Eq. (2) is unable 
to account for such effects. 
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D. NN interaction in the final state 



Effects of the NN interaction in the final state in various production reactions were investi- 
gated in detail in Refs. [33,34,51-53]. In particular, in Ref. [34] it was demonstrated that the 
NN FSI cannot be factorized from the production amplitude if one wants to obtain reliable 
quantitative predictions for the cross sections. This conclusion confirms the arguments given 
in Ref. [33] . Also it was shown in Ref. [34] that the use of the Jost function of some realistic 
NN potentials for the evaluation of FSI effects is invalid and may lead to results consid- 
erably different from the ones based on a proper calculation. In this work NN FSI effects 
are treated consistently, i.e. the calculation of loop integrals is performed with taking into 
account the full dynamics both from the NN T-matrix and from the production amplitude 
(cf. Eq. (1) or Fig. (2), respectively). In the calculation we utilize not only the CCF [47] 
NN model, i.e. the one which is used for the ISI as well, but also Bonn B [54] in order to 
examine the sensitivity of the results to the NN interaction model. Both models describe 
equally well the NN phase shifts in the 1 S'o and 3 Si partial waves, which are relevant for 
the r) production, for energies below the pion production threshold, i.e. in the energy region 
relevant for the FSI. 

Note that the Coulomb interaction in the pp final state is taking into account, applying the 
prescription described in Ref. [55]. 

III. RESULTS OF THE CALCULATION OF THE REACTION NN -► NNr] 

Our results for the reactions pp — > ppn, pn — > pni] and pn — > drj are presented in Fig. 8 and 
compared with empirical information from the Uppsala/Celsius [6,10-12] and Jiilich/COSY 
[8,9] accelerator facilities. The dashed curves correspond to the calculation based on the 
original MN model (cf. section II A) for the elementary 77-production amplitude and the 
CCF NN model [47] for the ISI and FSI. 

Evidently, the calculation based on the original Jiilich MN model [18] yields a qualitative 
overall description of the experimental data. This has to be certainly considered as success 
because, as mentioned before, there are no adjustible parameters in our calculation of NN — > 
NNr]. As one can observe from Fig. 8, for the ppr\ case we overestimate the cross section 
by approximately 30% whereas for pn — > pni] the model calculation underestimates the 
experimental data by about 50%. The situation for the reaction pn — > drj is very similar to 
pn — > pnrj since both processes are governed by the 1 = isospin channel. As is known from 
earlier investigations [26,32] for the latter reaction the contribution of the 1 = 1 channel is 
much smaller than the one for 1 = 0. 

However, it is interesting to investigate the sensitivity of our results to the MN amplitude 
and, specifically, to see whether the model predictions for NN — > NNr/ can be improved by 
modifications in the MN model via introducing some additional graphs, as described in sec- 
tion II A. For the original MN model the amplitudes of the ir- and p-exchange contributions 
to 77 production turned out to be almost orthogonal to each other. Thus, there is practically 
no interference between these contributions. The additional diagrams which we introduced 
into the MN model allow a slight modification of the orientation of these amplitudes in the 
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complex plane and, as a consequence, now interference effects do occur. Specifically, it is 
possible to generate a destructive interference in the isotriplet channel (1 = 1) and a con- 
structive interference in the 1 = case. This leads to a slight reduction of the cross section 
in the reaction pp — > ppn and to a significant enhancement for pn — > pnrj, cf. the solid lines 
in Figs. 8, bringing now the results quite close to the experiment. Also the prediction for 
the reaction pn — > dn is now in good agreement with the experimental data. 

A. Contributions of individual meson exchanges. 

Let us examine our model in more detail and analyze the contributions of the individual 
exchange diagrams. We do this for the calculation based on the extended MN model. The 
main differences between the results obtained for the original MN model and its extended 
version will be discussed at the end of this section. 

Corresponding results are presented in Fig. 9 (pp — > ppn) and in Fig. 10 (pn — > pnrj for 
1 = 0), where we show the production cross section for the individual meson exchanges 
together with the full results. Furthermore, in Table II we compare cross sections of the full 
calculation with results obtained without ISI, i.e. where in Eq. (1) was set to zero, at 
the specific excess energy Q = 35 MeV. The latter allows us not only to expose the relevance 
of the various meson exchange (rescattering) contributions for the n production cross section 
but also to elucidate that their relative importance is strongly influenced by the ISI. For 
example, in the model calculation without ISI the dominant contributions come from n and 
p exchanges, both being of comparable magnitude. However, the ISI reduces the p-exchange 
contribution much more strongly then the one from n exchange (cf. Table. II) leaving ir 
exchange as the only dominant production mechanism. The contributions from the other 
production mechanisms (a exchange, rj exchange, etc.) are already comparably small before 
including the ISI and they are also significantly reduced (though not as much as the p) by 
introducing the ISI. 

Let us emphasize in this context that phenomenological treatments like those based on 
Eq. (2) are unable to account for such effects resulting from details of the dynamics. They 
lead only to an overall reduction of the cross section independent of the production mech- 
anisms, as is indicated by the last 3 columns of Table II. (Note that the reduction factors 
used in Ref. [32] differ from those for CCF because the ones in the former work are based 
on the NN phase shifts listed at the SAID library [44]. The reduction factors of Ref. [26] 
are not obtained from Eq. (2) but from a different prescription.) 

We also want to mention that the ISI effects seen in our investigation are at variance with 
those reported in the only other study of 77-meson production where the ISI was taken 
into account explicitly [24]. In that work by Batinic et al. it was found that the combined 
distortions from the FSI and ISI are not sensitive to the dynamics of the production operator. 
We don't have an explanation readily at hand for this discrepancy. We can only conjecture 
that it might be due to differences in the employed MN — > nN transition amplitudes. 

As we have discussed above, after inclusion of the ISI ir exchange plays the dominant role in 
our calculation of the reaction pp — > ppn (see Figure 9). Nevertheless, the contributions of the 
other mesons, especially of the n and p exchanges and of the direct term are still significant 
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due to interference effects with the amplitude of the tt exchange. The n exchange and the 
direct term contribute constructively to the reaction pp — > ppn whereas the interferences 
from the a and p mesons are destructive. 

Note that the actual effect of those interferences depends also on whether the ISI is included 
or not. That can be best seen by the fact that the ISI reduces the individual meson-exchange 
contributions by factor 4 or more whereas the total cross section is only reduced by roughly 
a factor 2, cf. Table II. 

Let us now come to the reaction pn — > pni]. The cross section of this reaction is primarily 
determined by the isovector particles (n and p) in the 1 = channel, because their contri- 
butions are weighted by the large isospin coefficient (9 versus 1 for 1=1). In addition, for 
the relevant NN partial wave in the initial state, the 1 P\, it turned out that the variation 
of ISI effects for the various production mechanisms is not as pronounced as in the 1=1 
( 3 P ) case. Indeed, except for pion exchange, in this channel the reduction due to the ISI is 
almost universal as can be seen in Table II. 

The relative contributions of the individual meson exchanges to the reaction pn — > pnn (in 
the 1 = channel) are shown in Fig. 10 for the full calculation. As mentioned already, the 
dominant role belongs again to the pion exchange, though the contribution from p exchange 
is now much less suppressed as in the 1 = 1 case (cf. Fig. 9). The individual contributions 
from the other production mechanisms are again much smaller. Like in case of pp — > ppn 
interferences play a role once we add all contributions coherently. Specifically, the n- and the 
p-meson exchange exhibit opposite features as in the pp — > ppn case, i.e. the n is destructive 
while the p is constructive. The direct term acts again constructively and its contribution 
nearly cancels with the one resulting from n exchange. The contribution from a exchange 
turns out to be negligible. 

The interference effects between n and p meson exchange amplitudes for the original and 
extended MN models are illustrated in Table III, where we present the ratio of the coherent 
sum of the two amplitudes to the incoherent sum. For the original MN model this ratio is 
close to 1 (for both isospin channels) verifying that there are only small interference effects. 
For the extended MN model those interferences are much more pronounced. 

Table II lists also the phenomenological reduction factor that follows from Eq. (2) for the 
employed NN initial state interaction (in the column with the header 'CCF'). Thus, we can 
compare directly the reduction of the cross section that follows from the explicit inclusion 
of the ISI with the one suggested by the phenomenological prescription. It is evident that 
the prescription Eq. (2) doesn't work that well for the 7r and p exchange contributions. For 
the n and specifically for the a exchange, however, the results are fairly similar. 

B. Role of the r]NN coupling constant 

Now we would like to comment on the influence of the nNN coupling constant. As men- 
tioned before, in our calculation we take the value g^ NN /(ATr) = 1.8 that follows from SU(6) 
symmetry. This value is already much smaller than those used in OBE versions of the Bonn 
NN model [43,54], say, which are around 3 to 7. On the other hand, some studies of n 
production in NN collisions suggest that only still smaller nNN coupling constants allow 
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to describe the experimental data [26,30]. E.g., a value of only about g^NN/i.^ 71 ) = 0-4 was 
employed in [30] while in [26] the n contribution was even completely neglected. Therefore, 
it is interesting to see whether our results on n production based on the original MN model 
could be also improved by simply varying the nNN coupling constant in the production 
operator. 

In our model calculations g V NN enters via the n exchange contribution but also in the direct 
term. The latter is not negligible and interferes constructively with the dominant n exchange 
contribution in both isospin channels of the reaction NN — > NNn. On the other hand, as 
mentioned before, the n exchange contribution interferes constructively for the reaction 
pp — > ppn but destructively for the 1 = part of pn — > pni]. Indeed, the direct term and the 
77 exchange basically cancel each other in the dominant isospin 1 = channel of the reaction 
pn — > pni] and, therefore, the nNN coupling constant influences primarily the cross section 
of the pp — > ppn reaction. Thus a reduction of g^NN would only decrease the cross section 
for pp — > ppn. It would not alter the results for pn — > pnn and pn — > dn and, therefore, does 
not lead to an improvement for the latter two reaction channels. 

C. Sensitivity to the NN interaction. 

In order to examine the sensitivity of our results to differences in the NN interaction in 
the final state we also performed a calculation where the OBE potential Bonn B [54] was 
utilized. Corresponding results are presented in Fig. 11 by dashed-dotted lines. The results 
with the CCF iViV model are also shown (solid lines) in order to facilitate a comparison. As 
one can see from this figure, for the pp — > ppn channel the predictions of both NN models lie 
basically on top of each other. However, for the pn — > pnn channel the results based on the 
Bonn B potential are about 20% larger and, indeed, are practically in agreement with the 
experiment. At the same time the cross section for the pn — > dn channel is also enhanced 
by about 20% and now slightly overshoots the data. 

In any case, the differences between the results for the two considered AW models are not 
that large. But let us emphasize here an interesting by-product of the above comparison. 
Obviously it is not possible to achieve a simultaneous description of all three measured n 
production channels. Either we have agreement for the reactions pp — > ppn and pn — > dn 
(for CCF) and then pn — > pnn is off, or pp — > ppn and pn — > pnn are reproduced (Bonn 
B) and then the results for pn — > dn deviate from the data. Let us remind the reader that, 
close to threshold, there are only two independent amplitudes that determine those three 
reaction channels. The same incompatibility is also seen in the results of Faldt and Wilkin, 
cf. Fig. 5 of Ref. [26], and, therefore, we believe that it is not due to the specific production 
mechanism employed in our study. It remains to be seen whether this discrepancy will 
disappear when contributions from higher partial waves are included. In any case, we want 
to mention that a similar incompatibility exists also for near-threshold pion production, cf. 
the results in Fig. 3 of Ref. [56] and also in Ref. [57]. 
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D. Comparison with other model calculations. 



In this subsection we want to discuss qualitatively the differences between our model and 
other model calculations in the literature. We restrict the comparison to such calculations 
where also pp as well as pn induced n production channels are considered, i.e. to the works 
of Gedalin et al. [21] and Nakayama et al. [32], where results for the reactions pp — > pprj 
and pn — > pni] are presented, and to Ref. [26], where in addition the channel pn — > dn is 
considered. 

The relatively large value of 6.5 for the cross section ratio <J pn ^ pnr] / '<J PP ^ PPV established in 
the CELSIUS experiment [10] revealed that the dominant 77-production mechanism has to 
be of isovector nature. (Exchange of scalar mesons would yield a ratio of roughly one!) 
Furthermore, simple estimations [10,26], taking into account the different strengths in the 
final (pp versus pn) state interactions, strongly suggest that a production operator involving 
only a single isovector-meson exchange (-7T or p) might still fall short of describing the experi- 
ments and, thus, further production mechanisms should be relevant. Indeed, a scenario with 
a dominant contribution from isovector meson exchange is realized in all the models cited 
above, including ours. However, while p exchange is the dominant n production mechanism 
in Refs. [21,26] and in the alternative model of Ref. [32] it is tt exchange that provides the 
bulk of the cross section in our model and in the regular model of Nakayama et al. [32]. 
We should mention, however, that in our model 7r and p exchange are of pretty much the 
same strength and provide comparable cross sections before the ISI is included, as can be 
seen from Table II. If we would ignore ISI effects altogether, as done in Ref. [21], or use the 
phenomo logical prescription employed in Refs. [26,32] this feature would persist. 

The models differ also significantly when it comes to the role played by the smaller contribu- 
tions. For example in the model of Faldt and Wilkin the 7r exchange is very important. Its 
strong interference with the dominant p exchange contribution is the main mechanism which 
allows those authors to achieve a consistent description of the experimentally observed cross 
sections for the reactions pp — > ppn and pn — > pnn. On the other hand, in the alternative 
model of Nakayama [32] it is the interference of the u exchange with the dominant p ex- 
change that is responsible for obtaining a reproduction of the data. Finally, in the regular 
model of the latter publication agreement with the data is achieved by the interference of 
the 77 exchange with the dominant 7r exchange. 

In our model calculation based on the extended MN model it is the constructive interference 
between the ir- and p-meson contributions which yields the main enhancement in the cross 
sections for pn — > pnn and pn — > dn and brings the results close to the experiment. Thus, 
for the 1 = dominated reactions the mechanism is similar to the one in the model of Faldt 
and Wilkin. However, this is no longer true for the reaction pp — > ppn. Here, in our model 
the 7r-p interference does not play an important role. The reason for that are the strong 
effects from the ISI, discussed above, which reduce the p-exchange contributions much more 
strongly than those coming from ir exchange, cf. Table II and, therefore, suppress also the 
interference. The reproduction of the data for the reaction pp — * ppn is mainly due to a 
constructive interference of the dominant 7r-exchange contribution with the contributions 
from n exchange and from direct n production. Thus, for the 1 = 1 case the mechanism is 
similar to the one in the regular model of Nakayama [32] . 
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In the work of Gedalin et al. [21] several interference patterns are studied. But, in general, 
their results do not agree that well with the experiments in the energy region Q < 50 MeV 
considered. Because of those reasons we refrain from a more detailed comparision with 
their results. Let us mention, however, that the most prominent feature in their model 
is the absence of any interferences between the contributions of vector meson exchange 
and pseudoscalar (or scalar) meson exchange. The orthogonality of those contributions is 
achieved by a peculiar choice of interaction Lagrangians for the vector mesons. However, it 
seems to us that those Lagrangians do not fulfil the requirement of time reversal invariance. 
Also, they seem to be in contradiction to the ones applied by these authors in their own 
earlier investigations [20]. 

Interestingly, also in our calculation based on the original MN model there is practically no 
interference between the 7r and p exchange contributions, as already mentioned above. But 
in our model the origin of this feature is quite different. Here the orthogonality of the two 
amplitudes is generated by the dynamics of the coupled-channels nN model. In particular, 
the Born term pN — > iV*(1535) — > i]N that is the main source of the orthogonality in the 
model of Gedalin et al. is not even present in the original MN that we use in our study 
of i] meson production. Once this diagram is included in the coupled-channels MN model 
(cf. the extended MN model described in subsect. II A) interference effects become more 
pronounced rather than suppressed. 

IV. SUMMARY 

We performed a detailed theoretical calculation of the different channels of the reaction 
NN — > NNr] (pp — > ppr], pn — > pni], pn — > dr]) in the near-threshold region, i.e. for excess 
energies up to about 50 MeV. The production mechanisms which have been included consist 
of re-scattering terms with M = it, p,rj,a meson exchanges and the direct 77 production. 
Effects of the final and initial state interaction between the nucleons are fully taken into 
account. The calculation utilizes the CCF NN model for the treatment of the FSI and ISI 
and a realistic coupled-channels model of the tcN system for the evaluation of the MN — ■> 
r/N transition amplitudes. A qualitative agreement of calculated cross sections with the 
experimental data is achieved for all considered r\ production channels. This has to be 
certainly considered as success because there are no adjustable parameters in this model 
calculation. It is also shown that even a quantitative description of the data can be obtained 
if one introduces small modifications of the MN — > r/N amplitudes by exploiting some 
freedom in the pN — > r]N transition potential of the original MN model. 

In our model the dominant role in ^-meson production near threshold belongs to the re- 
scattering mechanism with intermediate pion exchange. The contributions from other meson 
exchanges, specifically from p and 77 as well as from the direct 77 production are smaller. 
However, these mechanisms are still important and have an influence on the cross sections 
due to their interference with the amplitude corresponding to the n exchange. 

Our study shows that ISI as well as FSI effects are not universal but depend, among oth- 
ers, on the concrete meson production mechanisms. Thus, a consistent treatment of these 
effects is very important for a quantitative comparison of model calculations with available 
experimental data on the considered reactions. Specifically, the interaction in the initial 
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NN system plays a crucial role since it affects the magnitude of the cross sections, leading 
to a significant reduction. 

The presented model should be viewed as a first step towards a consistent description of rj 
production in NN collisions and in meson induced reactions. Future work should avoid some 
of the technical approximations, e.g. by an explicit treatment of the three-body singularities 
[50] . In addition, and most importantly, further mechanisms for 77-meson production should 
be explored. As mentioned, 77-meson production via 00 rescattering is still missing in the 
present model but should be taken into account. Furthermore, rescattering contributions 
involving the A isobar in the two baryon intermediate states [56,58] should be investigated. 
Finally, higher partial waves should be included in order to facilitate an extension of the 
model calculation to higher excess energies and also to make a comparison with measured 
angular distributions, polarizations and invariant mass spectra meaningful. 
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TABLES 



TABLE I. Meson masses, coupling constants and cutoff masses utilized in the calculation. 
The values are taken from the cited references. Monopole type form factors are used at all 
mesons-baryon vertices. 





m M (MeV) 


2 

9mnn 


^mnn (MeV) 


Ref. 


7T 


138.03 


14.4 


1300 


[43] 


P 


769 


0.84 ( Kp = 6.1) 


1400 


[43] 


V 


547.45 


1.8 


1500 


[18] 


a 


550 


5.689 


1700 


[43] 



TABLE II. Contributions of individual meson exchanges to the AW — > NNrj cross section and 
the influence of the AW initial state interaction. The values correspond to the excess energy Q 
= 35 MeV. The column labelled 'ratio' exemplifies the actual reduction of the cross section after 
inclusion of the ISI. For comparison we also present phenomenological reduction factors employed 
in the model calculations of Nakayama et al. [32] and Faldt and Wilkin [26]. In the column 'CCF' 
we give the phenomenological reduction factor that follows from Eq. (2) for the employed AW 
initial state interaction. 







&NN—*NNr) [pb] 








phenomenological 
reduction factors 


pp — > pprj 


with ISI 


without ISI 


ratio 


CCF 


Ref. [32] 


Ref. [26] 


IT 


3.11 


9.35 


0.33 




T 


T 


T 


P 


0.25 


10.3 


0.02 










V 


0.24 


3.01 


0.08 




0.13 


0.19 


0.59 


a 


0.36 


3.17 


0.11 










full result 


5.69 


10.4 


0.55 




I 


I 


1 


pn — > pnr] (1=0) 


with ISI 


without ISI 


ratio 


CCF 


Ref. [32] 


Ref. [26] 


7T 


10.4 


26.12 


0.4 




T 


T 


T 


P 


4.30 


25.96 


0.17 










V 


0.16 


0.85 


0.19 




0.30 


0.27 


0.53 


a 


0.27 


1.05 


0.26 










full result 


24.48 


70.39 


0.35 




I 


I 


1 
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TABLE III. Illustration of interference effects between the reaction amplitudes based on the tt 
and p meson-exchange contributions. The values correspond to the excess energy Q = 35 MeV. 





Original model 


Extended model 




\M V +M p y 


\Mk+M p \ a 


\M^+\M P \^ 


\M^+\M P \* 


1=1 


0.88 


0.66 


1=0 


1.13 


1.45 
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FIG. 1. Production mechanisms for the reaction NN — > NNr] taken into account in our model: 
(a) rj production via MN — > rjN rescattering; (b) direct rj production. 
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FIG. 2. Diagramatic representation of our DWBA calculation. ^2^3 is the elementary 
^-production amplitude. The filled ellipses stand for the NN interaction in the final and initial 
states. 
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FIG. 3. The ttN phase shift and inelasticity 1 — rj 2 for the S\\ partial wave. The dashed curves 
represent the results of the original 7T./V model of Ref. [18,39] whereas the solid lines are the results 
of the extended model that includes the additional diagrams of Fig. 4. 
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FIG. 4. Additional diagrams containing the pNN* coupling which are included in the extended 
model. 
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FIG. 5. S-wave cross section for the reaction ir~p — > r/n. The dashed curve represents the 
results of the original ttN model of Refs. [18,39] whereas the solid line is the result of the extended 
model that includes the additional diagrams of Fig. 4. Data are taken from Ref. [41]. 



21 



50 




T lab (MeV) T lab (MeV) 



FIG. 6. Phase shift S and inelasticity parameter p (rj = cos 2 p) calculated for the 3 Pq partial 
wave using the CCF model [47]. The squares represent experimental phase shifts, extracted from 
the SAID library [44] and recalculated from the Arndt-Roper convention [48] to satisfy the condition 
for the S-matrix: S = rje 2zS . The arrows indicate the rj production threshold. 
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FIG. 7. Phase shift 5 and inelasticity parameter p ( 77 = cos 2 'p) calculated for the l P\ partial 
wave using the CCF model [47]. Same description of squares as in Fig. 6. 
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FIG. 8. Total cross sections of the reactions pp — > pprj (top), pn — > pnr? (middle) and — * drj 
(bottom). The dashed lines represent the results of our calculation with the original MN model 
of Ref. [18] , whereas the solid lines are based on the extended MN model described and discussed 
in Sect. II A. Data are taken from Refs. [6,8-12]. 
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FIG. 9. Contributions of the individual meson exchanges to the reaction pp — > pprj. The 
presented results are based on the extended MN model. The solid line is our full result, i.e. when 
all contributions are summed up coherently. 
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FIG. 10. Contributions of the individual meson exchanges to the reaction pn — > pnr/ (I = 0). 
The presented results are based on the extended MN model. The solid line is our full result, i.e. 
when all contributions are summed up coherently. 
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FIG. 11. Total cross sections of the reactions pp — > pprj, pn — > pnr/, and pn — > dr/ employing 
different NN models for the final state interaction. The solid lines represent the results with the 
CCF NN model [47] whereas the dashed-dotted lines were obtained for the Bonn B model [54]. 
The calculations are based on the extended MN model. 
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